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ABSTRACT 

We present a catalog of 5324 massive stars in the Small Magellanic Cloud (SMC), with accurate 
spectral types compiled from the literature, and a photometric catalog for a subset of 3654 of these 
stars, with the goal of exploring their infrared properties. The photometric catalog consists of stars 
with infrared counterparts in the Spitzer SAGE-SMC survey database, for which we present uniform 
photometry from 0.3 - 24 ^m in the UBVIJHKs+lRAC+MlPS24: bands. We compare the color 
magnitude diagrams and color-color diagrams to those of stars in the Large Magellanic Cloud (LMC), 
finding that the brightest infrared sources in the SMC are also the red supergiants, supergiant B[e] 
(sgB[e]) stars, luminous blue variables, and Wolf-Rayet stars, with the latter exhibiting less infrared 
excess, the red supergiants being less dusty and the sgB[e] stars being on average less luminous. 
Among the objects detected at 24 ^m in the SMC are a few very luminous hypergiants, 4 B-type 
stars with peculiar, flat spectral energy distributions, and all 3 known luminous blue variables. We 
detect a distinct Be star sequence, displaced to the red, and suggest a novel method of confirming Be 
star candidates photometrically. We find a higher fraction of Oe and Be stars among O and early-B 
stars in our SMC catalog, respectively, when compared to the LMC catalog, and that the SMC Be 
stars occur at higher luminosities. We estimate mass-loss rates for the red supergiants, confirming 
the correlation with luminosity even at the metallicity of the SMC. Finally, we confirm the new class 
of stars displaying composite A & F type spectra, the sgB[e] nature of 2dFS1804 and find the FO 
supergiant 2dFS3528 to be a candidate luminous blue variable with cold dust. 

Subject headings: catalogs- galaxies: individual (SMC)- infrared: stars- stars: early-type- stars: 
emission-line. Be- stars: massive 



1. INTRODUCTION 



^ Institute of Astronomy & Astrophysics, National Observa- 
tory of Athens, I. Metaxa & Vas. Pavlou St., P. Penteli, 15236 
Athens, Greece; bonanos@astro.noa.gr 

^ Space Telescope Science Institute, 3700 San Martin Drive, 
Baltimore, MD, 21218, USA; lennon, massa, sewilo, panagia, 
kgordon, meixner, shiao@stsci.edu 

^ European Space Agency, Research and Scientific Support 
Department, Baltimore, MD 21218, USA 

Department of Physics and Astronomy, Heidelberg Univer- 
sity, Albert-Ueberle-Str. 3-5, D-69120 Heidelberg, Germany; 
koehlinger@stud.uni-heidelberg.de 

^ Astrophysics Group, Lennard-Jones Laboratories, Keelc 
University, Staffordshire ST5 5BG, UK; jacco@astro.keele.ac.uk 
® UK Astronomy Technology Centre, Royal Observa- 
tory Edinburgh, Blackford Hill, Edinburgh, EH9 3HJ, UK; 
chris.evans@stfc.ac.uk 

IN AF/Osservatorio Astrofisico di Catania, Via Santa Sofia 
78, 1-95123 Catania, Italy 

Supernova Ltd., OYV #131, Northsound Road, Virgin 
Gorda, British Virgin Islands 

® Department of Astronomy, 475 North Char- 
ter Street, University of Wisconsin, Madison, WI 
53706, USA; brian@sal.wisc.edu, s_bracker@hotmail.com, 
meade@astro.wisc.edu 

Steward Observatory, University of Arizona, 933 North 
Cherry Avenue, Tucson, AZ 85721, USA; mblock, cengelbracht, 
kmisselt@as. arizona. edu 

Harvard-Smithsonian Center for Astrophysics, 60 Garden 
Street, MS 67, Cambridge, MA 02138, USA; jhora, tro- 
bitaille@cfa.harvard.edu 

Department of Astronomy, University of Virginia, PO Box 
3818, Charlottesville, VA 22903, USA; remy@virginia.edu 

Space Science Institute, 4750 Walnut Street, Suite 205, 
Boulder, CO 80301, USA; bwhitney@spacescience.org 



The Spitzer Space Telescope Legacy Surveys SAGE 
("Surveying the A gents of a Galaxy' s Evolution" 
Meixner et all [20061) and SAGE-SMC (IGordon et al.l 
2010f ) have for the first time made possible a com- 
parative study of the infrared properties of massive 
stars at a range of metallicities, by imaging both the 
Large and Small Mag ellanic Clouds (LMC and SMC). In 
B onanos et al.l (j2009l hereafter Paper I) , we presented in- 
frared properties of massive stars in the LMC (at 0.5 Z©). 
The motivation for that study was twofold: (a) to use the 
infrared excesses of massive stars to probe their winds, 
circumstellar gas and dust, and (b) to provide a tem- 
plate for studies of other, more distant, galaxies. Paper I 
was the first major compilation of accurate spectral types 
and multi-band photometry from 0.3—24 ^m for massive 
stars in any galaxy, increasing by an order of magnitude 
the number of massive stars for which mid-infrared pho- 
tometry was available. The recently completed SAGE- 
SMC survey offers the opportunity to extend the study 
of infrared properties of massiv e stars to a metallic ity of 
approximately 0.2 Zq (see e.g. iHunter et al.ll2007D . and 
enables the investigation of their dependence on metal- 
hcity, at least over the range 0.2-0.5 Zq. 

Infrared excess in hot massive stars is primarily due 
to free fre e emission fr om their ionized , line-driv en, stel- 
lar wi nds. IPanagia feTFellil (|1975j) and lWright fc Barlowl 
(|1975f ) first computed the free-free emission from ion- 
ized envelopes of hot massive stars, as a function of the 
mass-loss rate (M) and the terminal velocity of the wind 
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(woo)- The properties of niassive stars, and in particu- 
lar their stehar winds (which affect their evolution) are 
expected to d e pend on metallicity (Z). For example, 
iMokiem et al.l (|2007[ ) found empirically that mass-loss 
rates scale as M ~ ^o. 83±o.i6^ in good agreement with 
theoretical predictions (|Vink et al.ll2001[ ). The expecta- 
tion, therefore, is that the infrared excesses of OB stars 
in the SMC should be lower than in the LMC, given that 
M is lower in the SMCF^. Furthermore, there is strong 
evidence that the fraction of classical Be stars among B- 
type stars is higher at lower me tallicity (possibly due t o 
faster rotatio i i, as m easured bv lMartavan et aLll2007b[ ). 
IGrebel et al.l (fl99l were the first to find evidence for 
this, by showing that the cluster NGC 330 in the SMC 
has the largest fraction of Be stars of any known cluster 
in the G alaxy, LMC or SMC. M ore recent spectroscopic 
surveys (jMartavan et al.l [20101) have reinforced this re- 
sult. In Paper I, we showed that the Be stars in the LMC 
are easily discriminated by their mid-infrared colors; we 
are therefore interested in a comparison with the SMC 
to quantify the global dependence of the Be star fraction 
on metallicity. The incidence of Be/X-ray bina ries is also 
much higher in the SMC than in the LMC (|Liu et al.l 
I2005D , while the incidence of Wolf-Rayet (WR) stars is 
much lower; therefore, a comparison of infrared excesses 
for these objects is also of interest. Finally, there is in- 
terest in the metallicity dependence of the mid-infrared 
colors of red supergiants, which probe circumstellar dust. 

Following the same strategy as in Paper I, we com- 
piled a catalog of massive stars, which we cross-matched 
in the SAGE-SMC database to study their infrared prop- 
erties. This paper is organized as follows: Section [2] 
describes our spectroscopic and photometric catalogs of 
massive stars in the SMC, Section|3|presents the resulting 
color-magnitude and two-color diagrams, and Section 21 
presents the infrared excesses detected in various types of 
massive stars. Section [5] investigates the mass-loss rates 
in red supergiants, and Section|6|summarizes our results. 

2. CATALOG OF MASSIVE STARS IN THE SMC 

As in Paper I, we have compiled a catalog of massive 
stars with known spectral types in the SMC from the lit- 
erature. We then cross-matched the stars in the SAGE- 
SMC database, after incorporating optical and near- 
infrared photometry from recent surveys of the SMC. 
The resulting photometric catalog was used to study the 
infrared properties of the stars. In Section 12.11 we de- 
scribe the spectral-type catalog compiled from the liter- 
ature, in Section 12.21 the existing optical, near-infrared 
and mid- infrared surveys of the SMC, which were in- 
cluded in the SAGE-SMC database, and in Section lOl 
the cross-matching procedure and the resulting photo- 
metric catalog. 

2.1. Spectral-type Catalog 

The criteria for inclusion in the spectral type catalog 
of massive SMC stars were the availability of accurate 
coordinates and accurate (mainly optical) spectral clas- 
sifications that correspond to stars with masses > SM©. 
The 2dF Survev lEws et al.l (2dFS;l2004a|, 4161 objects) 

The dependence of i >oo on Z is negligib le for this range of Z 
(see ILeitherer etaniTQQ^ : lEvans et ani2004bl ) . 



comprises the largest catalog of SMC stars with spectro- 
scopic classifications and therefore makes up the back- 
bone of ou r spectral type catalosF^. Previous compila- 
tions (e.g. lAzzopardi et al.l llQTSt lAzzopardi fc Vigneaul 
[T979l [l982h were included with updated coordinates from 
Brian Skiff's listsj^ and spectral types from the liter- 
ature, when available. The literature search resulted 
in 5324 entries. We do not claim completeness; how- 
ever, not only have we targeted studies of hot massive 
stars (defined by their spectral types) in individual clus- 
ters and OB associations, but also studies of particular 
types of massive stars. The largest of t hese include: the 
survey of B and Be stars in NGC 330 (jMartavan et al.l 
I20?y7allbl . 314 stars), the VLT-FLAM ES Survey of Mas- 
sive Stars in NGC 330 and NGC 346 (jEvans et al.ll2006l 
232 stars), and unpublished additional FLAMES obser- 
vations in NGC 346 (PI: Evans; I. Hunter et al. 2010, 
in preparation, 244 OB & A- type stars). We have also 
added the 3 known luminous blue variables (LBVs), 5 su- 
pergiant B[e] stars (|Zickgrafll2006( ) . 69 Be/X-ray binaries 
(with spectral types from the updated online catalo£3 of 
uzova fc Popovll2005l l ■ early- type eclipsing an d spec- 
trosc opic binaries, an d red s uperg iants tMassev fc OlsenI 
I200I iLevesque et 311120061 mm 107 stars). The com- 
pleteness of the catalog depends on the spectral type, 
e.g., it is ~ 4% for the unevolved O stars in our catalog 
(out of an estimated total of 2800 unevolved stars with 
masses > 20 M^ in the SMC: Massc v 2009), while the 
WR catalog (jMassev et al.ll2003l ) is thought to be close 
to complete. 

Table [1] presents our catalog of 5324 massive stars, 
sorted by R.A., listing: the star name(s), coordinates in 
degrees (J2000.0), the reference and corresponding spec- 
tral classification. The names of the stars are taken from 
the corresponding reference. The spectral classifications 
in the catalog are typically accurate to one spectral type 
and one luminosity class. We retained about a dozen 
stars from the above studies with only approximate spec- 
tral types (e.g., "early B"). For double entries, we se- 
lected the most precise spectral classification available, 
usually corresponding to the most recent reference. 

2.2. Optical and Infrared Surveys of the SMC 

Several large optical and infrared photometric catalogs 
of the SMC have recently become available, enabling us 
to obtain accurate photometry for its massive star pop- 
ulation in the wavelength range 0.3 — 8 /im and in some 
cases up to 24 iim . The optical s urveys are as follows: the 
UBVR catalog of lMassevI ()2002[ ) with 85,000 bright stars 
in an area covering 7.2 square degrees, the U BVI Magel- 
lanic Clouds Photometric Survey (MCPS: IZaritskv et al.l 
I2OO2I ) including 5.1 million stars in the central 18 square 
degrees of the SMC, and the OGLE III catalog contain- 
ing VI p hotometry of 6.2 mil lion stars covering 14 square 
degrees (iUdal ski et al.ll2008f). The an gular resolution is 
- 2'.'6 for the catalog of lMassevI (f200^ . - 1"5 for MCPS, 
and - 1'.'2 for OGLE III. 

The existing near-infrared photometric catalogs in- 
clude: the Two Micron Ah Sky Survey (2MASS; 

By retaining all SMC stars from 2dFS, we are possibly includ- 
ing some stars with masses < SMq. 

ftp:/ /ftp. lowell.edu/pub/bas/starcats/ 
http://xray.sai.msu.ru/~raguzova/BcXcat/ 
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iSkrutskie et~al] 120061 extended by 6X2 MASS) and the 
targeted IRSF survey (jKato et al.ll2007D . which contains 
2.8 milhon sources in the central 11 square degrees of the 
SMC. 2MASS has a pixel scale of 2'.'0 pixel'^ an aver- 
age seeing of 2'.'5 and limiting magnitudes of J = 15.8, 
H = 15.1 and Ks = 14.3. IRSF has a pixel scale of 
0'.'45 pixel-i, average seeing of 1'.'3, 1'.'2, I'.'l in the JHKs 
bands, respectively, and limiting magnitudes of J = 18.8, 
H = 17.8 and Kg — 16.6. In the mid-infrared, the 
Spitzer SAGE-SMC survey uniformly imaged the whole 
SMC (--30 square degrees) in the I RAC and MIPS bands 
on two epochs in 2007-2008 (jGordon et al.ll2010D . The 
survey has recently produced a combined mosaic cata- 
log of 1.2 million sources. IRAC, with a pixel scale of 
1'.'2 pixel^^, yields an angular resolution of 1'.'7 — 2'.'0 
and MIPS at 24 /xm has a resolution of 6". 

Given the variation in the depth, resolution and spa- 
tial coverage of these surveys, we included the available 
photometry for the massive stars in our catalog from all 
the MCPS, OGLE III, 2MASS, IRSF and SA GE-SMC 
catalo gs. MCPS has incorporated the catalog of lMassevI 
(j2002l ) for bright stars common to both catalogs. Pho- 
tometry of higher accuracy, particularly in the optical, 
exists in the literature for many of the stars in our cata- 
log; however, it was not included in favor of uniformity. 

2.3. Photometric Catalog 
2.3.1. Matchtng Procedure 

We used the following SAGE-SMC data products to 
search for the mid-infrared counterparts to the 5324 mas- 
sive stars listed in Table [TJ the vl.O IRAC Catalog 
containing IRAC Epoch 1 data bandmerged with All- 
Sky 2MASS and 6X2MASS catalogs, and the MIPS 24 
/xm Epoch 1 Catalog. The details about the SAGE- 
SMC da ta processing and Le gacy data products can be 
found in lGordon etall (IMol) and in "The SAGE-SMC 
Data Description: Delivery 1" document available at the 
NASA/IPAC Infrared Science Archive Web siteEE. The 
IRAC Catalog and MIPS 24 ^m Catalog contain high- 
reliability sources and are subsets of more complete, but 
less reliable source lists: the IRAC Archive and MIPS 24 
/im full list (see the data delivery document for a com- 
parison between the catalogs and archives/full lists). 

We started constructing the photometric catalog for 
the massive stars in our list by selecting their SAGE- 
SMC IRAC Epoch 1 Catalog counterparts. We per- 
formed a conservative neighbor search with a 1" search 
radius, and selected the closest match for each source. 
We found mid-infrared counterparts for 3654 out of 
5324 sources. The IRAC Epoch 1 catalog (IRACC), 
MIPS 24 /im catalog, IRSF, MCPS, and OGLE III cata- 
logs were crossed-matched in the SAGE-SMC databascEB 
to provide photometry for sources over a wavelength 
range from 0.3-24 ^m. We used this "universal cata- 
log" to extract multi-wavelength photometry for IRAC 
sources matched to the massive stars. Specifically, 
for IRAC sources with one or more matches in other 
catalogs (all but 5), we only considered the closest 
matches between sources from any two available catalogs 
(IRAC-MIPS24, IRAC-MCPS, MIPS24-MCPS, IRAC- 
IRSF, MIPS24-IRSF, IRAC-OGLE III, MIPS24-OGLE 

1* http://irsa.ipac.caltech.edu/data/SPITZER/SAGE/ 
http:/ /mastweb. stsci.edu/scasjobs/ 



III), with distances between the matched sources of <1". 
For example, for a match between the IRAC, MIPS 24 
/xm and MCPS catalogs (IRAC-MIPS24-MCPS), we ap- 
plied these constraints on the IRAC-MIPS24, MIPS24- 
MCPS, and IRAC-MCPS matches. If the match was not 
the closest one or the distance was > 1", we dropped the 
pair. These stringent criteria were used to ensure that 
sources from individual catalogs for each multi-catalog 
match refer to the same star. Table [2] shows the break- 
down of the matched stars to the catalogs, such that 40 
stars were only matched to the IRACC, 214 only to the 
IRACC+IRSF catalogs etc. 

2.3.2. Catalog Description 

Table [3] presents our final matched catalog of 
3654 stars, with the star name, IRAC designation, 
UBVIJHKs+l^kC+MlVS2'i photometry and errors, 
reference paper, corresponding spectral classification and 
comments, sorted by increasing R.A.. Overall, the pho- 
tometry is presented in order of shortest to longest wave- 
length. The 17 columns of photometry are presented 
in the following order: UBVI from MCPS, VI from 
OGLE III, JHKs from 2MASS, JHR^ from IRSF, IRAC 
3.6, 4.5, 5.8, 8.0 ^tm (or [3.6], [4.5], [5.8], [8.0]) and MIPS 
24 jjLTn (or [24]). A column with the associated error 
follows each measurement, except for the VI photome- 
try from OGLE III. Henceforth, JHKs magnitudes refer 
to 2MASS photometry, whereas IRSF photometry is de- 
noted by a subscript, e.g. Jirsf- All magnitudes are 
calibrated relative to Vega (e.g. see iReach et al.l I2005L 
for IRAC). In Table |31 we summarize the characteristics 
of each filter: the effective wavelength Aoff , zero magni- 
tude flux (in Jy), angular resolution, and the number of 
detected stars in each filter. 

The spatial distribution of our 3654 matched sources 
is shown in Figure [l] overlaid onto the 8 /xm image of 
the SMC. For clarity, only a third of the early and late-B 
stars are shown; these are equally distributed along the 
SMC bar and wing. Other types are concentrated along 
the bar, partly due to the fact that most spectroscopic 
surveys have targeted the bar. The 8 /xm emission, which 
maps the surface density of the interstellar medium, is 
stronger in the SMC bar than the wing, since the former 
contains more star-forming regions. 

3. INFRARED COLOR-MAGNITUDE AND TWO-COLOR 
DIAGRAMS 

We divide the matched stars into 10 categories accord- 
ing to their spectral types: O stars, early (BO— B2.5) 
and late (B3— B9) B stars (the majority of these have 
supergiant or giant luminosity classifications), spectral 
type A, F and G type (AFG) supergiants, K and M red 
supergiants (RSGs), WR stars, supergiant B[e] (sgB[e]) 
stars, confirmed LBVs, Be/X-ray binaries and stars with 
composite features of both A and F stars (AFcomp; 
lEvans et all l2004al see also Section 0|. Note, that 
813 stars not included in any of the above categories are 
classified as "other" (and only shown in certain plots0. 
Most are A stars with luminosity class II or F/G stars 
withou t a luminosity classification from lEvans et al.l 
(|2004a|) . although with Vradiai > 100 km s"\ To allow 

20 Note that 5 Bextr, 5 Be, 1 Be(Fell), 1 Bpec, and 1 A7-F5e 
star are found among these. 
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for a direct comparison with the LMC resuhs of Paper 
I, we present the same color-magnitude and two-color 
diagrams (CMDs and TCDs), identifying stars in the 10 
categories described aboveFI 

Infrared [3.6] versus [3.6] - [4.5] and J- [3.6] CMDs for 
all the stars in the catalog are shown in Figures [2] and [3l 
The locations of all the SAGE-SMC catalog detections 
are represented by the gray two dimensional histogram 
(Hess diagram). The red giants form the clump at 
[3.6] > 15 mag, while the vertical blue extension contains 
late-type SMC and foreground stars (free-free emission 
causes the OB stars to have redder colors). The asymp- 
totic giant branch stars are located at [3.6] 11 mag, 
and J - [3.6] ~ 2 mag. The RSGs, sgB[e], and LBVs 
are among the brightest stars at 3.6 /im and occupy dis- 
tinct regions in the diagrams (as in the LMC). Most of 
the O and B stars are located along a vertical line at 
[3.6] - [4.5] ~ 0, as expected. The RSGs have "blue" 
colors because of the depression at [4.5] due to the CO 
band (see e.g. iVerhoelst et al1l2009f ). The sgB[e] stars 
are 0.6-0.8 mag redder and can have similar brightnesses 
to the RSGs. The late-B stars are brighter than the 
early-B stars because most of the former are luminous 
supergiants. A reddening vector for E{B — V) = 0.2 mag 
(correspondin g to some of the largest values for color ex- 
cess found bv lMassev et al.l [19951 : ILarsen et al1l2000[ ) is 
shown in Figure [3] to illustrate the small effect of red- 
dening, which decreases at longer wavelengths. We note 
that the reddest of the bright AFG s upergiants is the FOI 
eclipsing binary R47 (|Prieto et al.|[2008i ). In Figure [3l 
we find 3 early B-type stars with colors similar to the 
sgB[e] stars, although they are fainter at [3.6]. In order 
of increasing J — [3.6] color, these stars are: 2dFS0402, 
2dFS2411 and 2dFS2673 (all classified as B0-5V). Their 
evolutionary status is not clear, though they are worthy 
of future study (see discussion in Section l4?5]) . A compar- 
ison with the LMC CMDs reveals that sgB[e] stars are 
on average less luminous; while the RSGs are less dusty 
(see Section [SJ. 

In Figures |3] and [51 we present [8.0] and [24] versus 
[8.0] — [24] CMDs, respectively. At these wavelengths, 
only stars with strong [24] emission are detected, as the 
sensitivity of Spitzer drops sharply, while the flux from 
hot stellar photospheres also decreases. Stars without 
cool dust are located at [8.0] — [24] ~ 0, while dusty stars 
(including RSGs, LBVs, sgB[e] stars) are reddened. The 
faint population of sources with [8.0] — [24] ~ 3.5 mag 
corresponds to background galaxies. We detect dusty 
LBVs, sgB[e], RSGs, yellow supergiants and some dusty 
early-B stars in these CMDs. The detection of the fol- 
lowing stars - labeled as "other" - at [24] (in order of in- 
creasing magnitude at [24]): 2dFS352^3 (FO), 2dFS0712 
(G8), and 2dFS1829 (G2), provides further evidence for 
their supergiant luminosity class and thus membership 
to the SMC. The early B-type stars appearing in these 
diagrams, in order of increasing magnitude at [24] are 
as follows: 2dFS0402, 2dFS2673 (both of type B0-5V; 
discussed above), AzV9 (BOIII), AzV201 (Bl(Be-Fe)), 
AzV409 (or 2dFS2201; B0.5Ib), AzV230 (BOIb) and 

The conversion to absolute magnitudes in all CMDs is based 
on a true SMC distance modulus of 18.91 mag IHilditch et al.l 

[2003) . 

22 2dFS3528 is discussed further in Section [13^ 



AzV390 (B2Iab:). Of these, both AzV409 and AzV230 
have spectral energy distributions (SEDs) with disk-like 
characteristics in all Spitzer passbands (i.e. similar to the 
Be star SEDs shown in Paper I, Figure 17), although the 
former has a slight upturn at [24]. The nature of these 
objects is discussed in Section 14.51 Finally, the late-B 
and AFG supergiants present in these diagrams are hy- 
pergiants whose spectra and SEDs indicate the presence 
of substantial stellar winds. A comparison of the colors 
of various types of massive stars in the LMC versus SMC 
CMDs reveals that the colors of the RSGs differ the most 
(see Section |4. 5. 3p . 

In Figures El [71 and [51 we present TCDs using the near 
and mid-infrared photometry from our matched cata- 
log. We label stars according to their spectral types and 
overplot all the SAGE-SMC detections in gray as Hess 
diagrams. Stars without dust have mid-infrared colors 
near 0. Late-type stars can have J — Kg excesses of < 2 
mag, while the group of stars with colors K^ — [8.0] ~ 3, 
[4.5] - [8.0] - 2 and [8.0] - [24] ~ 3.5 mag correspond to 
young stellar objects and background galaxies. We have 
overplotted 6 simple theoretical models to guide the in- 
terpretation of the stars in these diagrams, which are 
described in detail in Paper I: (i) a blackbody (BB) at 
various temperatures, (ii) a power-law model cx v'^ , 
for —1.5 < a < 2, (iii) an OB star (represented by a 
30,000 K BB) plus an ionized wind, (iv) an OB star plus 
emission from an optically thin H II region, (v) an OB 
star plus 140 K dust, (vi) a 3,500 K BB plus 250 K 
dust. The most conspicuous stars in all TCDs are the 
sgB[e] stars, which have large excesses of ~ 4 mag in the 
K — [8.0] color, similar to the sgB[e] stars in the LMC. 
R4, a BO[e]-type LBV, has mid-infrared colors interme- 
diate between sgB[e] and Be stars. We find a remarkable 
bimodal distribution of the OB stars, most clearly seen in 
Figure [6l with the redder grouping corresponding to Oe 
and Be stars (see Section [^7T] for more details). Among 
these outliers are the Be/X-ray binaries and other emis- 
sion line stars, including some yellow supergiants. The 
near- infrared excesses found among the sgB[e] and WR 
stars are on average lower than in the LMC, presum- 
ably due to the lower metallicity of the SMC. Finally, 
the RSGs, in contrast to their distribution in the LMC 
TCDs, are found to be clustered, due to the lower dust 
production rates at the metallicity of the SMC (see Sec- 
tion [5]). 

4. INFRARED EXCESSES OF MASSIVE STARS IN THE 
SMC VERSUS THE LMC 

In this section we study the infrared excesses of specific 
spectral types of massive stars, comparing them with the 
infrared excesses found in Paper I for the same types of 
stars in the LMC. 

4.1. 0/Oe and early-B/Be stars 

In FigureslSlandfTUl we plot Jirsf versus Jirsf—[5-Q], 
JiRSF — [5.8] and Jirsf — [8.0] colors for the 208 
O and 1967 carly-B stars in our catalog, respectively, 
denoting their luminosity classes, binarity and emis- 
sion line classification properties by different symbols. 
We compare the observed colors with colors of plane- 
parallel non-LTE TLUSTY stellar atmosphere models 
(jLanz fc Hubenvl 120031 . 12007f) of appropriate metallicity 
and effective temperatures. For reference, reddening vec- 
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tors and TLUSTY models reddened by E{B - V) = 0.2 
mag are also shown. We clearly detect infrared excesses 
from free-free emission despite not having dereddened the 
stars, as in the LMC. At longer wavelengths, the excess 
is larger because the flux due to free-free emission for 
optically thin winds remains essentially constant with 
wavelength. Fewer stars are detected at longer wave- 
lengths because of the decreasing sensitivity of Spitzer 
and the overall decline of their SEDs. We find that the 
majority of early-B supergiants in the SMC exhibit lower 
infrared excesses, when compared to their counterparts 
in the LMC (see Paper I), due to their lower mass- loss 
rates, although certain exceptions exist and deserve fur- 
ther study. 

The CMDs allow us to study the frequency of Oe and 
Be stars, given the low foreground and internal red- 
dening for the SMC. Our SMC catalog contains 4 Oe 
stars among 208 O stars (see Figure IH]), of which one is 
bluer than the rest. There are 16 additional stars with 
JiRSF — [3.6] > 0.5 mag and Jirsf < 15 mag (includ- 
ing all luminosity classes), whose spectra appear normal 
(although the Ha spectral region in most cases was not 
observed). We refer to these as "photometric Oe" stars 
and attribute their infrared excesses to free-free emis- 
sion from a short-lived, possibly recurrent circumstellar 
region, whose Ha emission line was not detected dur- 
ing the spectroscopic observations either because the gas 
had dispersed or because the region was optically thick 
to Ha radiation or the observation spectral range just 
did not extend to Ha. Given the expectation of lower M 
at SMC metallicity, we argue that such a region is more 
likely to be a transient disk rather than a wind. Assum- 
ing these are all Oe stars, we find a 10 ± 2% fraction of 
Oe stars among the O stars in the SMC. The error in 
the fraction is dominated by small number statistics. In 
contrast, there are 4 Oe and 14 "photometric Oe" stars 
(with Jirsf — [3.6] > 0.5 mag and Jirsf < 14.5 mag) 
out of 354 O stars in the LMC (despite the higher M at 
LMC metallicity), which yields a 5 ± 1% fraction of Oe 
stars among O stars in the LMC. 

Turning to the early-B stars, the most striking fea- 
ture in Figure [10] is a distinct sequence displaced by 
^ 0.8 mag to the red. A large fraction of the stars 
falling on this redder sequence have Be star classifica- 
tions, although not all Be stars reside there. Given that 
the circumstellar gas disks responsible for the emission 
in Be stars are known to completely vanish and reap- 
pear be tween spectra taken eve n 1 year apart (see re- 
view bvl Porter &: Riviniu£ll2003l and references therein), 
the double sequence reported here provides further evi- 
dence for the transient nature of the Be phenomenon. A 
bimodal distribution at the L— band was previou sly sug- 
gested by the study of iDoughertv et a"Lr()1994D . which 
included a sample of 144 Galactic Be stars. Our larger 
Be sample, which is essentially unaffected by reddening, 
and the inclusion of all early-B stars, clearly confirms 
the bimodal distribution. It is due to the much larger 
number of Be stars classified in the SMC, in comparison 
to the LMC, as well as the higher fraction of Be stars 
among early-B stars in the SMC, which is 19 ±1% versus 
4 ± 1% in the LMC when considering only the spectro- 
scopically c onfirmed Be stars (cf. ~ 17% f or < 10 Myr 
BO-5 stars: IWisniewski fc BiorkmanI [2006D . Excluding 



the targeted sample of iMartavan et all ()2007al lb[) does 
not significantly bias the statistics, since the fraction only 
decreases to 15 ± 1%. We caution that incompleteness in 
our catalogs could also affect the determined fractions, 
if our sample turns out not to be representative of the 
whole population of OB stars. 

We proceed to define "photometric Be" stars as early-B 
type stars with an intrinsic color J/hsf — [3.6] > 0.5 mag, 
given that a circumstellar disk or envelope is required to 
explain such large excesses. Including these "photomet- 
ric Be" stars and using the same color and magnitude 
cuts as for the "photometric Oe" stars above, yields frac- 
tions of Be stars among early-B stars of 27 ± 2% for the 
SMC and 16 ± 2% for the LMC (cf. 32% fr om young 
SMC clusters; IWisniewski fc BiorkmanI |2006[ ). Table E] 
summarizes these results. W e compare our result s with 
the fractions determined by 'Maede r et al.l ([T999) from 
young clusters, i.e. 39% for the SMC and 23% for the 
LMC, finding ours to be lower, although the sample se- 
lections were very different. 

These preliminary statistics (available for the first time 
for Oe stars) indicate that both Oe and Be stars are twice 
as common in the SMC than in the LMC. We emphasize 
the importance of including the "photometric Be" stars, 
which significantly increase the frequencies of Oe/O and 
Be/ear ly-B stars determined and are crucial when com- 
paring such stars in different galaxies. This novel method 
of confirming Oe and Be star candidates from their in- 
frared colors or a combination of the ir optical and in- 
frared colors, as recently suggested bv llta et"an (j2010aD 
is com plementary to the detailed spectroscopic analyses 
by e.g. lNegueruela et al.l (|2004f ) on individual Oe stars to 
understand their nature, although it is limited to galax- 
ies with low internal reddening. We finally note that the 
spectral types of Oe stars in the SMC (07.5Ve, 07Ve, 
04-7Ve and 09-BOIH-Ve) and the LMC (09Ve (Fe II), 
07:Ve, 08-9IIIne, 03e) are earlier than those of known 
Galactic Oe stars, which are all found in the 09-BO range 
(jNegueruela et al.ll2004D . 

To illustrate the double sequence more clearly, we 
present the same information in a histogram in Figure ITT] 
for the early-B stars in both the LMC (from Paper I) 
and SMC, divided into magnitude bins. For stars with 
Jirsf — [3.6] ~ 0, we find the brighter, more luminous 
stars to be redder (cf. top 3 panels of Figure ITT]) , as in 
Paper I. In the second and third panels, the mean values 
of the Jirsf — [3.6] colors for the (blue, red) peaks are: 
(-0.07, 0.74) and (-0.09, 0.71) mag for the SMC and 
(-0.01, 0.66) and (-0.06, 0.79) mag for the LMC. The 
blue peak for the SMC stars is thus bluer than for the 
LMC, which is primarily due to the higher foreground 
reddening toward the LMC. In the second panel we find 
the separation between peaks to be larger for the SMC 
(0.81 mag) than the LMC (0.67 mag), implying that the 
SMC Be stars have larger infrared excesses than their 
LMC counterparts. Unfortunately, in the third panel, 
the red peak of the LMC stars is barely discernible, there- 
fore we cannot confirm this difference. We also ran a 
Kolmogorov-Smirnov test to check whether the stars are 
drawn from the same population (assuming the sample 
is unbiased). This yielded a probability of 5. Ox 10"'^ (for 
the stars in the 2nd row) and 5.7x 10~^ (3rd row), imply- 
ing that the underlying populations are indeed different. 
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Finally, we note that the brightest Be stars in the 
SMC {JiRSF ^ 13.2 mag) are brighter than the bright- 
est Be stars in the LMC {Jirsf ~ 13.4 mag), i.e. 
there is a 0.7 mag difference in absolute magnitude, 
given the 0.5 mag diffe r ence in the distance moduli. 
IGarmanv fc Humphrevsl (jl985f) were the first to make 
such a comparison, finding the Be stars in the Magel- 
lanic Clouds to be up to 1 mag brighter when compared 
to Galactic Be stars, although their observation could 
be due to a selection effect, related to the selection of 
field versus cluster Be stars. Even though our sample is 
not dominated by such selection effects, we also find the 
brightest Be stars in the SMC to be intrinsically more lu- 
minous than their LMC counterparts. Given that winds 
are weaker at low metallicity, it is possible that there 
is a mechanism that allows disks to form at higher lu- 
minositiej^. Such a possibility would have interesting 
implications when extrapolating t o even lower metallici- 
ties (see, e.g. lEkstrom et al.ll2008l ). 

4.2. Be/X-ray Binaries 

iDravl ()2006[ ) found that both metallicity and a star 
formation burst must be invoked to explain the large 
number of high mass X-ray binaries in the SMC, and 
consequently the large number of Be/X-ray binaries. Of 
the 69 Be/X-ray binaries in our SMC photometric cat- 
alog, 21 (or 30 ± 8%) were matched in the SAGE-SMC 
database (versus 4 matches out of 20 Be/X-ray binaries 
or 20±11% in the LMC). Interestingly, some SMC Be/X- 
ray binaries are brighter than their LMC counterparts, 
despite being farther away, in agreement with the find- 
ing that Be stars in the SMC occur at higher luminosity. 
A comparison with the infrared colors of the Be/X-ray 
binaries in the LMC does not reveal any differences. The 
SEDs of the SMC Be/X-ray binaries show excess due to 
the disks around the Be stars, as in the LMC. Finally, the 
multiple VI and JHKg measurements reveal variability 
for many systems, which deserves further study. 

4.3. Wolf-Rayet stars 

Of the 12 WR stars in our spectral-type catalog, 10 
were matched in the SAGE-SMC database. HD5980 
(SMC-WR5;AzV229 WN6h;LBV binary) is labelled as 
an LBV in all diagrams, and is the only WR star with a 
detection at [24] . We compare the colors of the WR stars 
to those of their LMC counterparts (see Paper I) and find 
that their infrared excesses are much lower, possibly in- 
dicating that their winds are much weaker than those of 
LMC WR stars. Specifically, we find the SMC WR stars 
to have "bluer" colors: —0.2 < J — [3.6] < 0.5 (versus 
0.2 - 1.7 for the LMC sample), -0.2 < J - < 0.1 
(versus - 0.6), 0.4 < Ks - [8.0] < 1.2 mag (versus 
0.4 — 1.8). The origin of the hydrogen absorption lines 
seen in 10 of the 12 WR stars (classified as "WN-|-abs" or 
"WNha" ) remains unclear, since radial ve locity surveys 
have only confirmed 4 as binary systems (jFoellmi et al.l 
|2003[) . The absorption in the other stars could alterna- 
tively be ex plained as photosph eric absorption (due to 
weak winds: lMassev et al.ll2003[) . The infrared photom- 
etry made available for the WR stars in this paper, will 

This would also explain the larger infrared excesses found in 
Be stars in the SMC, described in the previous paragraph. 



enable modeling of their SEDs to determine the proper- 
ties of their stellar winds and perhaps the origin of the 
absorption lines. 

4.4. AF composite stars 

lEvans et al.l (j2004aD discovered a class of objects 
with peculiar composite spectra (AF composites or AF- 
comp), which yield an A- type classification from their 
Ca II K lines, but an F-type classification from other 
metal lines and the G-band. For example the classi- 
fication AFA3kF8 denotes such an AF composite star, 
with Ca II K like A3 and G-band like F8. Our 
photometric catalog contains 16 out of the 20 AF 
compos ite stars from the 2dF Survey of lEvans et al.l 
(l2004al) . The AF composite stars that additionally ex- 
hibit forbidden emission lines (2dFS1804;AFA3kF0/B[e] 
and 2dFS2837;AFA5kF0/B[e]) are labeled and hence- 
forth treated as sgB[e] stars. 

The infrared colors of the AF composite stars provide 
further clues to their nature and confirm their astrophys- 
ical peculiarity. In Figure [21 their [3.6] — [4.5] colors are 
indistinguishable from AFG supergiantj^. while in Fig- 
ure [3] they clearly have red colors, with 0.5 < J — [3.6] < 
1 mag. None are detected at [24] . In Figure IH all 4 AF 
composite stars with [8.0] detections show near-infrared 
colors in the range 0.3 < J — Kg < 0.9 mag; thus, in- 
termediate between blue and red supergiants, possibly 
indicating these stars are on blue loops or in interact- 
ing binaries. The SEDs indeed are quite "flat" ; thus, not 
consistent with normal single stars, while the multiple VI 
and JHKs measurements reveal variability for several 
sources, which in some cases can be confirmed b y light 
curves available from th e OGLE database ( Udalski et al.l 
IT997t ISzvmafiskil l2005l) . We note that 2dFS2945 (HV 
11519) is a known 6 Cep pulsator. 

4.5. Dusty stars 

Our catalog includes 166 stars with detections at [24], 
which indicate the presence of dust. Of these, 97 are 
RSGs, 5 sgB[e] and 3 LBV stars, all exhibiting dust, as 
expected. Luminous supergiants with strong winds or 
disk-like SEDs, also exhibit infrared excess at [24] and 
[8.0] and therefore appear in Figures |3] and O There are 
4 BO-5 V type stars: 2dFS0402, 2dFS2411, 2dFS2673, 
and 2dFS3701 that have very fiat SEDs. Additionally, 
we find a class of 44 dusty OB stars (16 late-0 and 28 
early-B stars, making up 2% of our sample), with sim- 
ilar characteristics to the peculiar p o pulat ion of dusty 
OB stars reported by iBolatto et all ()2007l 190 09-B3 
stars, corresponding to 5% of their sample) in their S ^MC 
survey and recently confirmed bv llta et al.l (|2010bf ) with 
data from AKARI. Hereafter, we refer to these stars as 
"dusty OB stars". They are characterized by a photo- 
spheric SED out to [5.8] and strong [24] excess, which 
cannot be attributed to free-free emission. Most have 
main-sequence luminosity classes, although there are a 
few giants and supergiants. In Figure [T^l we plot SEDs 
for the 4 stars with unusually flat or rising SEDs, 2 repre- 
sentative examples from among the 44 "dusty OB stars" 
(2dFS1172, 2dFS1192) and 2 examples of luminous su- 
pergiants with infrared excess due to strong winds or 

24 The brightest AF composite star at [3.6] is 2dFS5049. 
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disks (AzV230, AzV409;2dFS2201). For completeness, 
in Figure [131 we show 2 examples of late-B supergiants 
with [24] detections and the 5 AFG supergiants that have 
[24] detections. Note, the FO supe rgiants R47 (a rare 
eclipsing binary; iPrieto et aL '20 08) and 2dFS3528 (an 
LBV candidate, see Section l4.5.2|) arc of particular inter- 
est. 

The dust detected in the class of "dusty OB stars" 
could either be directly related to the OB stars or sim- 
ply be contamination given the volume of the SMC sam- 
pled with the 6" resolution of MIPS. While inspection 
of their environments shows that they are located in 
or near young associations, they are frequently not the 
youngest objects (earliest types), e.g. NGC330 has 3 
"dusty OB stars" n earby and a turn-off of around B2 
(jLennon et al.lll993l ). We examined the spectra for 18 of 
these sources (12 from the 2dF Sur vey and 6 from the 
VLT-FLAMES Survev: lErons et al.|[20061 which aU ex- 
hibit sharp nebular-like Balmer emission lines; for those 
with red-optical spectra, 7 also have Ha nebular emis- 
sion and 3 have weak [N II] and [S II] lines, indicating the 
presence of a faint H II region. Many sources for which 
we do not have spectra available are either listed in Ha 
emission line catalogs or lie in regions of faint nebulosity. 
Given that most "dusty OB stars" are not found within 
very young regions, we suggest that these sources repre- 
sent unrelated cirrus hot spots or cases where hot stars 
are ionizing nearby dusty molecular clouds (rather than 
young stellar objects) producing both the nebular emis- 
sion lines and emission from warm dust (as in Cepheus 
B, see e.g. iPanagia fc fbl^lTQSl . 

4.5.1. Supergiant B[e] stars 

In the SMC photometric catalog, we have detected 
7 luminous sources with colors typical of sgB[e] stars 
(see Paper I for an introduction), i.e. M3 g < —8, 
[3.6] - [4.5] > 0.7, and J - [3.6] > 2 mag (see Figures [l 
and 131) . Five of these are previously known sgB [e] stars 
(with R50; B2-3[e] being the brightest in ah IRAC and 
MIPS bands), while R4 (AzV16) is classified as an LBV 
with a sgB[e] spectral type. In addition to these, we 
find that 2dFS1804 (AFA3kF0/B[e]) has a very similar 
SED (and therefore infrared colors) to the know n sgB[e] 
2dFS2837 (AFA5kF0/B[e]). lEvans et all (l2004ah also re- 
marked on the similarity of their spectra. We therefore 
confirm the supergiant nature of 2dFS1804 and present 
its SED in Figure [TH along with those of other sgB[e] 
and LBV stars. We note that both are photometrically 
variable. The similarity of the SEDs of these sgB[e] stars, 
despite the various optical spectral classifications, implies 
that all are the same class of object. The cooler, com- 
posite spectral types indicate a lower mass and perhaps 
a transitional stage to or from the sgB[e] phenomenon. 
The only difference we find between the sgB[e] stars in 
the SMC versus the LMC is that on average they are 
~l-2 mag fainter (in absolute terms). 

4.5.2. Luminous Blue Variables 

All 3 known LBVs in the SMC: R4 (AzV16, 
BO[e]LBV), R40 (AzV415, A2Ia: LBV) and HD5980 
(WN6h;LBV binary), were all detected at infrared wave- 
lengths. In all CMDs and TCDs, R4 is the more red- 
dened LBV. It occupies a location similar to other sgB[e] 
stars on the CMDs; however, in Figures |6l and [71 it lies 



between the sgB[e] and the Be stars, whereas the colors 
of HD 5980 (a well kno wn eccentric eclipsing binary, see 
e.g. iFoellmi et al1[2008l) are similar to those of the LBVs 
in the LMC. We plot their SEDs in Figure (Hand find 
them to differ, given their very different spectral types. 
Moreover, we find evidence for variability, which can be 
confirmed from existing light curves in t he All Sky Auto- 
m ated Survey (AS AS: lPoimahski|[2002[ ). as pointed out 
bv iSzczvgiel et all (|2010[ ). who studied the variability of 
the massive stars presented in Paper I in the LMC. The 
various SED shapes and spectral types observed depend 
on the time since the last outburst event and the amount 
of dust formed. 

Finally, we note that the FO star, 2dFS3528 (Skl66) 
also lies among the sgB [e] stars in Figures |4l |5l and [51 
and near R4 (AzV16; BO[e]LBV) in the other CMDs 
and TCDs, because it has a similar SED (cf. Upland Fig- 
ures [T31) . Although Sanduleak originally classified this 
star as "OB:", its 2dF spectrum is clearly F-type, but 
with strong Ha emission h aving broad wings, which indi- 
cate a wind. lOsmerl (|1973[ ) designated it as a nonmember 
based on its optical colors. It is included in the catalog of 
Masscy (2002]^ and indeed has some colors that are too 
blue, while others are too red for its s pectral type. How- 
ever, t he radial velocity measured by lEvans fc HowarthI 
(|2008l +U7 ± 8 km s^^) firmly places it in the SMC. 
The SED of 2dFS3528 is distinctly sgB[e]-like and most 
similar to that of R4. All the above characteristics, along 
with the variability detected in the near- infrared, make 
this star an LBV candidate. Further monitoring is highly 
desirable to confirm its LBV nature. 

4.5.3. Red Supergiants 

The RSGs are among the brightest mid-infrared 
sources in the SMC, with absolute magnitudes reach- 
ing — 12 mag at [3.6], and —13 mag at [8.0] and [24]. 
However, in the LMC, some RSGs occur at higher [24] 
luminosities. SMC_018592 (K0-2I) is the brightest RSG 
at [3.6], [8.0] and [24] (see Figures (H HO . In Figure p 
the following 3 RSGs have "bluer" colors than the rest (in 
order of increasing [3.6] magnitude): SMC_030135 (K2I), 
SMC_035231 (K2I) and SMC-049990; 2dFS1517. The 
latter has two cl assifications: K5I b y iMassev fc OlsenI 
(|2003D and G2 bv lEvans et"all (|2004aD . making it a vari- 
able candida te (given also its radial velocity: -1-143 ± 
5 km s-\ lEvans fc HowarthI I2008D . In Figure H 
SMC_035231 (K2I) is the outlier at (0.4, 0.6), making 
it a possible foreground dwarf; while SMC_058738 (K2I) 
is the outlier at (1.0, 0.5). We note that SMC_055188 
(M3-M4I) and SMC.046662 (MOI) were reported by 
ILevesque et al.l (|2007[) to have variable spectral types; 
they are included in our catalog and are among the more 
reddened RSGs. 

A comparison of the SMC CMDs and TCDs to the cor- 
responding ones of the LMC from Paper I, reveals that 
the distribution of RSG colors is quite different. The 
SMC RSGs have a dominant population of non-dusty 
RSGs and a population of dusty RSGs, with a distinct 
color distribution. In contrast, almost all RSGs in the 
LMC are dusty to some degree, with a rather continuous 
distribution in color. This difference is due to the lower 

25 [M2002] 80339; V = 12.92, B - V = 0.15, U ~ B = -0.45, 
and V - R = 0.17 mag. 
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dust content (cx Z), which is not compensated for by the 
fower wind velocity (oc -s/Z, see Ivan Lom]|20001 I2OO6I: 
Ivan Loon et al.ir20Q8[). Furth ermore, since both the dust 
production rate (cf. ISloan et al. 2008) and the dust-to- 
gas ratio depend on Z (|van Loonll2006[) . M is indepen- 
dent of Z, though th e wind velocities of R SGs dep end on 
metallicity (see e.g. iHabing et al.l 119941 Marshall et al.l 
12004^). While the optical and near- infrared colors of 
RS Gs have long been known to display a Z dependence 
(cf. lArpI 119591 : lOestreicher et al.llT997f l. the systematics 
of their mid-infrared col ors have only been investigated 
for limited samples (e.g. iGroenewegen et al.l l2"009). 

5. RED SUPERGIANT MASS-LOSS RATES 

Our large samples of RSGs in the SMC and LMC lend 
themselves to a study of the Z dependence of the mass- 
loss rate. Figure [15] presents the Mk^ versus Kg ~ [24] 
CMD for both SMC and LMC RSGs (from Paper I), 
on which the RSGs are spread out as a result of var- 
ious amounts of circumstellar dust. The fainter RSGs 
in the SMC concentrate around a vertical sequence at 
— [24] sa 0.35 mag, but all other RSGs have signifi- 
cantly redder Kg — [24] colors indicative of the emission 
from dust. A series of models for a star with spectra l type 
MO and luminosity 85,000 Lq (|G roenewegen 2006) are 
plotted for reference, showing that the dust remains op- 
tically thin at near-infrared wavelengths for — [24] < 5 
mag. 

The Kg ~ [24] color is a good measure of the opti- 
cal depth of the circumstellar dust envelope. We use 
the results from radi ative transfer modeling presented in 
iGroenewegenl (|2006l ) for a star of luminosity L — 3000 
Lq, wind speed = 10 km s~^, dust-to-gas ratio -0 ~ 
0.005, surrounded by silicate grains with a size of 0.1 /im 
and a condensation temperature of 1000 K. The relation 
between this reference M and Kg — [24] color is shown in 
Figure [HI for a star of spectral type MO (typical for the 
SMC stars) and again for a star of spectral type M6 (at 
the cool extreme of the LMC stars in our sample). For 
colors Kg — [24] < 0.8 mag, the relation is rather steep; 
given the inaccuracies in photometry and stellar parame- 
ters we will not assign mass-loss rates to these non-dusty 
stars. For colors 0.8 < Kg — [24] < 5 mag, the relation is 
quite well approximated by the formula: 

log(Afrcfcrcncc) - -8.6 + 0.5 X {Kg - [24]). (1) 

The deviations from this formula are not (much) more 
than a factor of two for individual stars, and will tend to 
average out when considering the sample as a whole. 

The true M is obtained by applying scaling relations 
valid for dust-driven winds (cf. Ivan Loonl[2007f ): 

M = M„ference(«/10)(0.005/?/;)(L/3000)°-5, (2) 

where the wind speed scales as (e.g. lMarshall et a"ni2004D 

V ^ 10(VV0.005)°•^(L/3000)°■2^ (3) 

For the dust-to-gas ratio ip, we adopt a value roughly in 
proportion to the metallicity of the stellar population, 
namely, ij ) k, 0.001 for the SMC and V ~ 0.002 for the 
LMC (cf. Ivan Loon|[2000l I2006D . The luminosity is ob- 
tained by scaling the band magnitude of the model 
star from iGroenewegen (,2006.) : for the SMC stars we 



adopt Kg — 8.1 mag for the 3000 Lq star, at a refer- 
ence distance of 8.5 kpc (this corresponds to a spectral 
type MO) and Kg = 7.9 mag for the LMC stars (spectral 
type about M2). The distance modulus to the SMC is 
taken as 18.91 mag, and to the LMC as 18.41 mag (to 
be consistent with Paper I). The inaccuracies in the lu- 
minosities resulting from this approximation are of the 
order of 20%, which is small for our purposes. 

The estimated mass-loss rates show a clear positive 
correlation with luminosity (see Figure 1171) , confirming 
previous results (^ir Loo n ct al. 1999, 2005; van Looij 
120061 : IGroenewegen et aD (20091 ) with a larger sample of 
stars. Apparently, the fainter RSGs are not very profi- 
cient in their mass loss. However, these mass-loss rates 
were derived by assuming that dust condenses in these 
winds at the efficiencies expected from their metal con- 
tent. This may not be true, particularly in the relatively 
warm atmospheres of early-M or late-K type stars, where 
the dust-to-gas ratio can b e much less than assumed here 
(cf. Ivan Loon et~aIll2005D . Indeed, at L < 9 x 10^ Lq 
most of the RSGs in the SMC show very little evidence 
of dust at all, but they may still exhibit a purely gaseous 
wind. On the other hand, the values for the wind speed 
are estimated to be 9-12 km s^^ for the SMC sample 
(12-19 km s~^ for the LMC sample). This is consis- 
tent with expectations for a fiuid in which the grains are 
well coupled to the gas and thus can transfer momen- 
tum from the radia tion field onto the bulk matter (cf. 
iMarshall et al.l[2004l ) ; much l ower values co uld bring this 
assumption into question (cf. iMcDonald et al. 2009 ). 

The mass-loss rate is next compared to the rate at 
which matter is converted into energy in the central en- 
gine of the star, to assess whether much mass is lost be- 
fore the star is taken elsewhere on its evolutionary path. 
In the context of core-helium burning RSGs, the most 
relevant process is the triple-a p rocess, which has an ef- 
ficien cy ?7 « 5.9 X 10^'' erg g~^ (jKippenhahn fc Weigerti 
[1990) :' we also compare to the CNO cycle, which governs 
hydrogen shell burning and has rj k, 6.1 x 10"'^* erg g"'^. 
The nuclear burning rate (in M© yr~^) is then: 

Afnuclcar = 6.087 X 10^(L/?7). (4) 

We find that essentially all RSGs in our samples lose 
mass through stellar winds at a significantly lower rate 
than that at which helium is burned in their cores (Fig- 
ure [TT]), often by an order of magnitude difference. This 
is in stark contrast to samples of RSGs that were se- 
lected from among th e brightest mid-infrared sources; 
Ivan Loon et al.l (|1999l ) showed that most of those RSGs 
have mass-loss rates similar to the nuclear burning rate, 
with a few extreme objects exceeding that rate by a sig- 
nificant amount. The optically selected RSGs in the 
SMC and LMC are therefore found to be in a more sta- 
ble phase of RSG evolution, in contrast to the RSGs 
selected as bright mid-infrared sources, which have a 
shorter lifetime. The prescriptions of Ide Jager et al.l 
(11988') for 3000, 4000 and 5000K coincide with the most 
luminous LMC RSGs, but are too high for the less lumi- 
n ous RSGs (in b o th ga laxies). The 3000K prescription 
of Ide Jager et al.l (|1988[ ) better reproduces the luminous 
RSGs in the SMC, whereas the warmer prescriptions bet- 
ter reproduce the luminous RSGs in the LMC. This is 
in conflict with the observed differences in temperature. 
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with the SMC RSGs having earher spectral types than 
their LMC counterparts. 

The SMC sample on the whole exhibits lower dusty 
mass- loss rates than the LMC sample (Figure [T7| , and 
RSGs in the SMC appear to contribute relatively more 
at higher luminosities (L > 9 x lO'* L©), while RSGs in 
the LMC also contribute noticeably at fainter luminosi- 
ties. Both observations are likely related to the warmer 
atmospheres (earlier spectral types) of the SMC RSGs 
compared to those in the LMC. We thus conjecture that 
dusty mass loss plays a less important role in the SMC 
than in the LMC, while other forms of mass loss may 
instead be more important. For instance, as recently 
suggested (|van Loon et al.|[2010aD . the warmer winds of 
metal-poor RSGs could be more susceptible to the cou- 
pling of Alfven waves, which could provide a driving 
mechanism in the absence of radiation pressure. 

We emphasize that our samples miss the dustiest 
RSGs, which are in brief episodes of very high mass loss; 
therefore, our results only apply to the optically bright 
RSG phase. A few such very dusty RSGs (known to 
exist in the Magellanic Clouds, e.g. 'Wood et al.' 1992'; 
Groenewegen fc Blomni acrt 1998; van Loon ct al. 1999, 
2010all bl: iSrinivasan et al.l 12009:) contribute more dust 
than the entire samples of optically bright RSGs con- 
sidered here. 

6. SUMMARY 

This paper presents the first catalogs of accurate spec- 
tral types and multi-wavelength photometry of massive 
stars in the SMC, which are used to study their infrared 
properties. The spectroscopic catalog contains 5324 mas- 
sive stars, with accurate positions and spectral types 
compiled from the literature, and constitutes the largest 
such catalog currently available for any galaxy. The pho- 
tometric catalog comprises uniform 0.3-24 /im photom- 
etry in the UBVIJHKs+lRAC+MlPS24 bands for a 
subset of 3654 stars that were matched in the SAGE- 
SMC database. The low foreground reddening toward 
the SMC, and the identical distance of the stars mini- 
mize systematic errors due to reddening and enable the 
investigation of infrared excesses. As in Paper I, we con- 
struct CMDs and TCDs, and discuss the position of O, 
early and late-B stars, WR, LBV, sgB[e], classical Be 
stars, RSGs, AFG supergiants, and Be/X-ray binaries 
on them. These diagrams are useful for interpreting in- 
frared photometry of resolved massive stars in nearby 
galaxies at low metallicity. 

A comparison of the infrared colors of the massive stars 
in the SMC to those of their counterparts in the LMC 
(presented in Paper I) reveals differences that are due 
to the different evolution at SMC metallicity. The main 
results of our study concern the emission line Oe and Be 
stars, and the RSGs. We detected a clear bimodal distri- 
bution of early-B stars, with the redder sequence corre- 
sponding to classical Be stars and therefore propose that 
Be stars (and similarly Oe stars) can be easily discrimi- 
nated photometrically, by their infrared colors. We find 
the fraction of emission line stars in the SMC (10 ± 2% 
for Oe and 27 ± 2% for Be) to be double that of the LMC 
(5 ± 1% for Oe and 16 ± 2% for Be), when including the 
"photometric" Oe and Be stars. This is the first time the 
frequency of Oe stars is determined beyond the Galaxy 
and at subsolar metallicities. We also find Be stars and 



Be/X-ray binaries to occur at higher luminosity, sgB[e] 
stars to be on average less luminous than their coun- 
terparts in the LMC (at [3.6]), and WR stars to have 
smaller excess, all due to the different evolution at the 
lower metallicity of the SMC. 

The infrared colors of the RSGs in the SMC are found 
in most cases to be consistent with little dust, with only 
the most luminous sources showing excess emission pre- 
sumably from circumstellar dust. This is in contrast 
to the generally dusty RSGs in the LMC, and agrees 
with the expectation that the dust content in metal- 
poor RSGs is lower. We find that the mass-loss rates 
in SMC RSGs correlate positively with luminosity, as in 
the LMC. Finally, we confirm the astrophysical peculiar- 
ity of the c omposite A & F type spectra discovered by 
lEvans et al . (2004a), the sgB[e] nature of 2dFS1804, and 
find the FO supergiant 2dFS3528 to be an LBV candi- 
date. This paper thus demonstrates the wealth of infor- 
mation contained in the SAGE-SMC survey, and enables 
studies of the infrared properties of massive stars as a 
function of metallicity, in combination with the SAGE 
survey of the LMC. Studies of particular classes of mas- 
sive stars and whole massive star populations in nearby 
galaxies, as outlined in Paper I, are obvious directions 
for follow-up. 
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Institute of Technology, funded by the National Aeronau- 
tics and Space Administration and the National Science 
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RA(J2000) 



Fig. 1. — Spatial distribution of massive stars with IRAC counterparts, overplotted onto the 8 fim SAGE image (Right Ascension and 
DecHnation are given in degrees). Different symbols denote different spectral types. For clarity, only a third of the early and late B stars 
are shown; these are equally distributed along the SMC bar and wing. Other types are concentrated along the bar, due to the fact that 
most spectroscopic surveys have targeted the bar. 
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Fig. 2. — [3.6] versus [3.6] — [4.5] color magnitude diagram for 
massive stars with IRAC counterparts in the SAGE database. The 
conversion to absolute m agnitudes is based on a true SMC distance 
modulus of 18.91 mag ijHilditch et al.l [20051 1. Different symbols 
denote different spectral types. The locations of all the SAGE 
detections are shown in gray as a Hess diagram. The RSGs, sgB[e] 
and LBVs are among the most luminous stars at 3.6 /im. 
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Fig. 4. — [8.0] versus [8.0] — [24] color magnitude diagram, with 
SAGE detections shown as gray points. The RSGs, sgB[e] and 
LBVs are also among the most luminous stars at 8 fim. The names 
and classifications of stars labeled as "other" are (in order of in- 
creasing [8.0] magnitude): 2dFS0712 {G8), 2dFS3528 (FO), and 
2dFS1829 (G2). 
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Fig. 3. — Same as Figure[2l but for the [3.6] versus J — [3.6] color 
magnitude diagram. The reddening vector for E{B — V) = 0.2 
mag is shown. The longer baseline separates the populations more 
clearly. 
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[8.0]-[24] 

Fig. 5. — Same as Figure g] but for the [24] versus [8.0] — [24] 
color magnitude diagram. The brightest sgB[e], RSGs and LBVs 
are among the most luminous stars at 24 fira. The names and 
classifications of stars labeled as "other" are (in order of increasing 
[24] magnitude): 2dFS3528 (FO), 2dFS0712 (G8), and 2dFS1829 
(G2). 
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Fig. 6. — J — Ks versus Ks — [8.0] diagram for massive stars in our catalog. The locations of all the SAGE detections are shown in gray 
as a Hess diagram. The solid lines represent models (described in Section|3ll: (i) a BB at various temperatures, as labeled, (ii) a power law 
model ¥i, oc i/", for —1.5 < a < 2, (iii) an OB star plus an ionized wind, (iv) an OB star plus emission from an optically thin H II region, 
(v) an OB star plus 140 K dust, (vi) 3,500 K blackbody plus 250 K dust (dashed line). The sgB[e], RSGs and the clump of "red" early-B 
stars that correspond to Be stars occupy distinct regions on this diagram. 
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K5]-[8.0] 

Fig. 7. — Same as Figure [6] but for the [3.6] — [4.5] versus [4.5] — [8.0] diagram. The majority of hot massive stars lie between the 
blackbody and OB star +wind model, illustrating that a BB is a good approximation in the infrared. 
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Fig. 8. — Same as Figure|6] but for the [3.6] — [8.0] versus [8.0] — [24] diagram. The more uniform distribution of RSGs, as compared to 
the LMC, is due to the lower dust content. 
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Fig. 9. — Infrared excesses (Jirsf versus JinsF — [3-6], Jirsf — [5-8] and JmsF — [8-0]) for 208 O stars. Supergiants are shown in 
yellow, giants in green, main-sequence stars in blue, stars with uncertain classifications ("other") in red, binaries with a large circle and Oe 
stars with an X. The solid lines correspond to 30kK and 50kK TLUSTY models with logg = 4.0. A reddening vector for E{B — V) = 0.2 
mag is shown, as well as reddened TLUSTY models by this same amount (dotted lines). The more luminous stars exhibit larger infrared 
excesses, which increase with A. 
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Fig. 10.— Same as Figure[9l but for 1967 early-B stars. The solid lines correspond to 20kK, logg = 3.0 and 30kK, logg = 4.0 TLUSTY 
models. A reddening vector for E{B — V) = 0.2 mag is shown, as well as reddened TLUSTY models by this same amount (dotted lines). 
The distinct redder sequence, corresponding to the Be stars, implies the possibility of identifying or confirming Be stars photometrically. 
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Fig. 11. — Histograms of the 3irsf — [3.6] colors of 1967 early B 
stars in the SMC (left) and 586 in the LMC (from Paper I, right), 
divided into magnitude ranges. The main peaks of the histograms 
shift to redder colors for brighter stars. The redder secondary peaks 
correspond to the Be stars, which are more numerous in the SMC. 
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Fig. 12. — SEDs of OB stars with [24] detections, normalized 
by their J— band fluxes (dashed line) and offset for display pur- 
poses. Normalized TLUSTY model atmospheres at SMC metal- 
licity (30kK, logg = 3.00 for the main sequence stars; 20kK, 
logg = 2.25 for the 2 supergiants) are overplotted as solid lines for 
comparison. The MCPS, IRSF and SAGE measurements (filled 
circles) and the 2MASS and OGLE measurements (open circles) 
are also connected by a solid line. The dotted curves correspond to 
TLUSTY models reddened by E{B-V) = 0.25 and 0.50 mag. Sev- 
eral classes of object are represented (see Section 14.511 . 2dF'S1172 
and 2dFS1192 are examples of dusty OB stars. Variability is 
present in several sources; however, in the case of 2dFS0402 it 
is likely due to a mismatch with a close infrared source. 
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Fig. 13. — SEDs of 2 late-B supergiants and of the 5 AFG super- 
giants detected at [24] , normalized by their J— band flux es (dashed 
line) and offset for display purposes, as in Figure [121 The SED 
of the F supergiant 2dFS3528 most resembles that of R4 (AzV 16, 
BO[e]LBV); we suggest it to be an LBV candidate. 
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Fig. 14.— SEDs of the LBVs and sgB[e], all detected at [24], 
normalized by their J— band flux es (dashed line) and offset for 
display purposes, as in Figure 1121 TLUSTY model atmospheres 
are overplotted for comparison; 15kK, logg = 1.75 for the B8 and 
later types; 20kK, logg = 2.25 for the earlier types. The last two 
SEDs with composite spectral types closely resemble the sgB[e] 
SEDs; note, these stars vary. 
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Fig. 15.— Mk, versus - [24] CMD of our sample of RSGs in 
the SMC (filled circles; this work) and LMC (open circles; Paper 
I). The fainter SMC RSGs cluster around Ks - [24] = 0.35 (dotted 
line). A model from IGroenewegen' (2000) for a star with spectral 
type MO and luminosity 85,000 Lq with increasing amounts of dust 
toward redder color is shown for reference. Most LMC stars and 
many SMC stars show clear excess emission due to dust, which is 
most pronounced at high luminosity. 
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Fig. 17. — M versus luminosity L for our samples of RSGs in the 
SMC (filled circles; this work) and LMC (open circles; Paper I). 
No mass-loss rates could be determined for stars with Ks — [24] < 
0.8 mag; these are shown in the bottom panel at their respective 
luminosities (crosses). Shown for reference are the nuclear burning 
rates for the triple-a process at work in helium burning in the cores 
of RSGs (dashed line), as well as the CNO c ycle in hydrogen shell - 
burning (dotted line) and the prescription of lde Jager et all II1988I) 
for 3000, 4000 and 5000K (solid lines, labeled by the first digit of 
the temperature). 




Fig. 16. — M as a function of Ks — [24] for two sets of models from 
IGroenewegen ( 2006) with spectral types MO (filled circles) and M6 
(open circles). The solid line indicates the approximation described 
by the labeled equation, which is valid for a star of luminosity 3,000 
L0, wind speed 10 km s~^ and dust-to-gas ratio 0.005, for the 
range 0.8 < Ks — [24] < 5 mag. We use scaling relations to apply 
it to the RSGs in the Magellanic Clouds. 
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TABLE 1 

Catalog of Spectral Types for 5324 SMC Massive Stars 



Star 


RA (J2000) 


Dec (J2000) 


Reference'' 


Classification 


Name" 


(deg) 


(deg) 




& Comments 


2dFS0001 


6.944625 


-73.27911 


E04 


09.5: 


2dFS0002 


7.036375 


-73.49611 


E04 


B0-5V 


2dFS0003 


7.089875 


-73.17531 


E04 


B0-5V 


2dFS0004 


7.102458 


-73.88650 


E04 


B0-5III 


2dFS0005 


7.137542 


-73.30739 


E04 


B0-5V 


2dFS0006 


7.243000 


-73.46847 


E04 


A2II 


2dFS0007 


7.246125 


-73.28139 


E04 


B1-3II 


2dFS0008 


7.322250 


-73.14731 


E04 


B0-5IV 



Note. — Table[T]is available in its entirety in the electronic version 
of the Journal. A portion is shown here for guidance regarding its 
form and content. 

Star designations: 'Evans et al.' {2dFS; '2004a'), 'Azzopa rdi et al.l 
(AzV ; 1975); Azzopardi & Vigncau (AzV; 1979, 1982), Feast ct al.l (R: 
[l960h , ,Sanduleah (Sk; 1968, 1969), Masscv ct al. (MFC; 1989). 
fa Ref erence k ey: ( A09) Antoniou ct al. (2009), (B95) Bar ba et af] 
(IT995I) . fgOl) iBartzakos ct al.. (.2001. ). (F03) .Foellmi c t al. f200a), 
( E04) [Ej ans ct al. (2004a), (E06) lEvans et al.l |2 006), (ElO) 

Slunter et al. (in Drcoaration), (GC87) 'Garmanv et al.ni987) , (H83) 
lumphrcys (1983), (H03) Harries ct al. (2003), (H05) Hildi tch et al.l 
2005), (L97) Lcnnon (1997 ), (L06) Lcvcsquc ct al. (2006), (L07) 
.cvcsQuc ct al. (2003), (MPG89) Masscv ct al. (1989), (M95) 
[Masscv ct al. (1995), (MOO) Masscv ct al. (2000), '(M02) Massp 
12002), (MO03) Masscv & Olscn (2003), (M03) Masscv ct al. (200^ 
(M05) Masscv ct al. (2005), (M06) Mcnnickcnt ct al. (2006), (M07) 
Martavan ct al. (2 007al lhl, (M09) Masscy ct al. (2009), (P08) 
Pricto ct al. (2008 ), (R05 ) Raguzova & Fopov (2005), (SB97) 
Smit h Ncubig & Br uhweileri (IT997I) . (W07) iWisniewski et al.l 8200711 . 
(Z06) Zickeraf {200(§'. 



TABLE 2 

Statistics for the 3654 Matched Stars 



Catalogs Matched # Stars 



IRACC 40 

IRACC+IRSF 214 

IRACC+IRSF+OGLE 352 

IRACC+MCFS 3 

IRACC+MCFS+IRSF 57 

IRACC+MCFS+IRSF+OGLE 2733 

IRACC+MCFS+OGLE 54 

IRACC+MCFS+MIPS24 3 

IRACC+MIFS24 6 

IRACC+MIFS24+IRSF 35 

IRACC+MIFS24+IRSF+OGLE 10 

IRACC+MIFS24+MCFS+IRSF 80 

IRACC+MIFS24+MCFS+IRSF+OGLE 32 

IRACC+OGLE 35 
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TABLE 3 

0.3-24 fim Photometry of 3654 Massive Stars in the SMC 



Star Name" 


IRAC Designation 


RA(J2000) 


Dcc(J2000) 


U 


f^(7 


B 




. Rof.'' 


Classification & Comments 


2dFS0005 


J002832.82-731827.0 


7.137542 


-73.30739 


16.33 


0.042 


17.063 


0.023 . 


. E04 


B0-5V 


2dFS0006 


J002858.15-732807.0 


7.243 


-73.46848 


16.648 


0.048 


16.605 


0.056 . 


. E04 


A2II 


2dFS0010 


J002923.69-734424.6 


7.349375 


-73.73997 


16.706 


0.041 


17.318 


0.033 . 


. E04 


B0-5V 


2dFS0016 


J003033.08-735709.9 


7.6385 


-73.95272 


15.977 


0.052 


16.051 


0.03 


. E04 


AOII 


2dFS0019 


J003049. 36-740000.3 


7.706292 


-73.99995 


16.006 


0.05 


16 


0.041 . 


. E04 


AOII 


2dFS0024 


J003107.03-732709.9 


7.780042 


-73.45267 


16.611 


0.114 


17.119 


0.034 . 


. E04 


B3IV 


2dFS0026 


J003111. 52-735936.1 


7.798458 


-73.99333 


15.744 


0.049 


16.292 


0.033 . 


. E04 


B1-3V 



Note. — Table [3l is available in its entirety in the electronic version of the Journal. A portion is shown here for guidance regarding its form 
and content. 

Star desi gnat ions: l Evans et al.l (2dFS:l2004al '), IAz2opardi et al.l (AzV: [T975I ): lAzzopardi fc Vigneaiil (AzV: [T979I . [rgsg) . IFeast et afl ('R: ll96a) , 
ISanduIeaU (Sk; 119681. 1 1969). Masscv et al.l (MPG: 119891) . 

^ Reference key: (B 95) Barba ct al. fl^ ), (BO ll [Bartzakos et al.l 120011'), (E OA) lEvans eFal] |[200i a |') (E06) lEvans et all 12003) , (ElO ) 
[Hjintcr et al. (iri preparation), (F03) FocUmi ct al. (2003), (GC87) Garmany ct al. (T987'). (H83)'Humphrcvs ( 19831. (H031 IHarries et al.l 12003') , 
(H05) Hilditch ct al. (2005), (L97) Lcnnon (1997), (L06) Lcvcsquc ct al. (2006), (L07) Lcvcsquc ct al. (2007), (MPG89'| iMassev et al] 09891'), 
(M95) Masscv ct al. (1995), (M02) Masscv (2002), (MO03) Masscv & Olscn (2003), (M06) Mcnnickcnt ct al. (20 03), (M07) | Martayan ct alj 
i2007a b), (M09) Masscy ct al. (2009), (P08) Pricto ct al. (2008), (R05) Rag uzova &: Pop ov (2005). (SB97) Smith Neubig fc Bruhwcilcr ( 1997. ), 
(W07) iWisniewski et al.i m07t] , (Z06) iZickgrai (|200Q1 . 



TABLE 4 

Filter & Detection Characteristics 





Filter 




Zero mag 


Resolution 


# stars 






(/im) 


flux (Jy) 


(") 


detected 


u 


0.36 


1790 


1.5/2.6 


2921 


B 


0.44 


4063 


1.5/2.6 


2962 


V 


0.555 


3636 


1.5/2.6 


2962 




I 


0.79 


2416 


1.5/2.6 


2844 


VoGLE 


0.555 


3636 


1.2 


3152 


lOGLE 


0.79 


2416 


1.2 


3166 


J 


1.235 


1594 


2.5 


3039 


H 


1.662 


1024 


2.5 


3007 


Ks 


2.159 


666.7 


2.5 


2789 


JlRSF 


1.235 


1594 


1.3 


3449 


HiRSF 


1.662 


1024 


1.3 


3413 


KsjRSF 


2.159 


666.7 


1.3 


3345 




3.6] 


3.55 


280.9 


1.7 


3509 




4.5] 


4.493 


179.7 


1.7 


3111 




5.8] 


5.731 


115.0 


1.9 


1038 




8.0] 


7.872 


64.13 


2 


509 




24] 


23.68 


7.14 


6 


166 



TABLE 5 
Oe and Be Star Fractions 



Fraction 



Oe/O 
(Oc+"phot Oc")/0 
Bc/early-B 
(Be+"phot Be")/early-B 



LMC SMC 



1 ± 0.6% 2 ± 1% 

5 ± 1% 10 ± 2% 

4 ± 1% 19 ± 1% 

16 ± 2% 27 ± 2% 



